Engineering superpositions of displaced number states of a trapped ion 
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We present a protocol that permits the generation of a subtle superposition with 2 e+1 displaced 
number states on a circle in phase space as target state for the center-of-mass motion of a trapped 
ion. Through a sequence of £ cycles involving the application of laser pulses and no-fluorescence 
measurements, explicit expressions for the total duration of laser pulses employed in the sequence 
and probability of getting the ion in the upper electronic state during the I cycles are obtained and 
analyzed in detail. Furthermore, assuming that the effective relaxation process of a trapped ion can 
be described in the framework of the standard master equation for the damped harmonic oscillator, 
we investigate the degradation of the quantum interference effects inherent to superpositions via 
Wigner function. 



I. INTRODUCTION 



In quantum mechanics, the nomenclature 'nonclassical states' has been employed by theoretical physics for approxi- 
mately seventy five years in order to designate the states of quantum systems whose statistical properties present 
genuine quantum effects (e.g., photon antibunching, sub-Poissonian statistics, squeezing, and quantum interference 
effects inherent to superposition states) without having analogous effects in classical mechanics [I]. During this period, 
different experimental techniques have been developed for generating and detecting both trapped and travelling 
nonclassical states, and more recently the nonclassical electronic and vibrational states of trapped ions. For instance, 
the density matrices and Wigner functions associated to the Fock states, thermal states, coherent states, squeezed 
vacuum states, and Schrodinger cat states (entangled position and spin superposition states), were reconstructed in 
laboratory through a beautiful experiment involving the quantum states of motion (motional states) of a harmonically 
bound 9 Be + ion [2-5]. After this experiment, a considerable number of papers dedicated to generation and detection 
of nonclassical motional states have appeared in the literature [6-19]. In particular, Matos Filho and Vogel [6] have 
considered a class of nonlinear coherent states (NCS) which exhibits interesting nonclassical features such as strong 
squeezing and self-splitting with pronounced quantum interference effects, and showed that they may appear as 
stationary states (also recognized by the authors as dark states) of the center-of-mass (CM) motion of a trapped and 
bichromatically laser-driven ion far from the Lamb-Dicke regime. Man'ko et al. [12] have extended the results obtained 
by Matos Filho and Vogel for NCS on a circle, where the influence of nonlinear effects on the Wigner functions was 
discussed in detail. Furthermore, Kis et al. [14] have introduced a method by which any pure state of the quantum 
harmonic oscillator can be represented in a limiting sense as a NCS, and showed through a physical example how 
to prepare a highly excited Fock state in an ion trap based on the concept of NCS. On the other hand, Moya-Cessa 
et al. [10] have shown how an arbitrary superposition of coherent states can be created on a line in phase space for 
the motion of a single trapped ion. Pursuing this line, Duarte Jose and Mizrahi [II] have proposed three schemes 
to engineer circular states (superposition of N coherent states on a circle in phase space) for the CM motion of a 
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trapped ion, whore the total duration of laser pulses and the probability of getting the ion in the upper electronic 
state were determined for each process. In addition, the authors also have shown how the interference effects between 
the components of the subtle superposition can produce the Fock states. 

Recently, many authors have investigated different sources of decoherence in experiments involving trapped ions 
and predicted interesting results which permit us to give reasonable explanations on the phenomenological decay 
rate of Rabi oscillations [20-23]. In this sense, Schneider and Milburn [20] have considered as decoherence source the 
intensity and phase fluctuations in the exciting laser pulses, and showed that a simple master equation description can 
be obtained since the stochastic processes involved are white noise processes. Serra et al. [21] have examined, through 
the analogy with the physics of surface electrons in liquid helium, the mechanism of damping and heating of trapped 
ions associated with the polarization of the residual background gas induced by the oscillating ions themselves. In 
particular, the authors have demonstrated that the decay of Rabi oscillations observed in experiments on 9 Be + can be 
attributed to the polarization phenomena. Budini et al. [22, 23] have assumed that the origin of decoherence of the 
nonclassical motional states is due to the coupling of the vibronic modes with classical fields and to the finite lifetime 
of the electronic levels, and showed that these interactions lead to a dispersive-like decoherence dynamics. On the 
other hand, Poyatos et al. [24] have shown how to design different couplings (due to the absorption of a laser photon 
and subsequent spontaneous emission) between a single ion trapped in a harmonic potential and environment. In this 
scheme, the variation of the laser frequencies and intensities allows one to 'engineer' the coupling and select the master 
equation which describes the motion of the ion. Turchette et al. [25] also have presented results from an experimental 
study of the decoherence and decay of motional states due the interaction with several types of engineered reservoirs. 
Now, independently of the damping mechanisms to be considered, the degradation of the quantum interference effects 
in superpositions of motional states is always verified. 

One of the most important applications using trapped ions was established by Cirac and Zoller [26] in the context 
of quantum computation, where the authors have shown that a set of N cold ions interacting with laser light and 
moving in a linear trap provide a realistic physical system to implement a quantum computer. The main features of 
this proposal are that (i) decoherence can be made negligible during the computation process, (ii) the implementation 
of n-bit quantum gates between any set of ions is relatively straightforward, and (iii) the quantum bit readout can be 
performed with efficiency approximately equal to one (a quantum bit or qubit refers to a two-state system characterized 
by {|0), |1)}). However, the experimental realization of a quantum computer requires isolated quantum systems acting 
as the qubits, and the presence of controlled unitary interactions between the qubits allowing the construction of 
the controlled-NOT gate (basically, a controlled-NOT is defined by the operation |ei) |e 2 ) — > |ei)|ei © £2) with © 
denoting addition modulo 2, and €1,2 = 0, 1). Thus, if the qubits are not sufficiently isolated from the environment, 
different mechanisms of decoherence can destroy the quantum interferences that make the computation. The first 
experimental implementation of a fundamental quantum logic gate that operates on prepared quantum states in 
experiments involving trapped ions was realized by Monroe et al. [27]. Following the scheme proposed in [26], the 
authors have demonstrated a controlled-NOT gate on a pair of qubits which illustrates the basic operations necessary, 
and the problems associated, in the construction process of a large scale quantum computer (in this experiment, the 
switching speed of the controlled-NOT gate is approximately 20 kHz and the decoherence rate is of a few kHz) . After 
the original ion-trap proposal of Cirac and Zoller, a number of modifications and extensions to their idea have appeared 
in the literature (e.g., see Refs. [28-33]). In summary, the investigation of noise sources in such promising quantum 
systems turns out to be a crucial step toward the implementation of a quantum logic processor, and consequently, of 
a quantum computer. 

According to Kis et al. [14]: "Nonclassical states of the electromagnetic field and the atomic center-of-mass motion 
have played an important role in recent years, due to their relation with fundamental problems in quantum mechanics 
and to the many possible applications, ranging from high-resolution spectroscopy to low-noise communication and 
quantum computation. However, the generation of these states is usually a demanding experimental challenge." In 
the present contribution, we propose a systematic scheme which permits us to engineer superpositions of displaced 
number states on a circle in phase space as target states for the CM motion of a trapped ion. These superpositions 
were studied by Marchiolli et al. [34] , where the authors have shown that (i) the interference effects among the state 
components present an analogy with diffraction patterns arising in an N slit Young- type experiment, and (ii) the 
interference and correlation effects are connected with the nondiagonal term of the quasiprobability distributions. In 
general, the superpositions of N displaced number states on a circle in phase space can be defined as follows: 



N N 

|*W(/3)> ee NW£D(ft)|n) = NW £ \n,(3 r ) (n e N) , 

r— 1 r— 1 



(i) 
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where D(/3 r ) = exp (/? r a^ — /3*a) is the displacement operator with (3 r — \/3\e Wr and 9 r+1 — 9 r = 2<p for <p € [0, 27r], 

{JV-l ~| _1 / 2 

N + 2j2 r e~ m2 sin2[(A, ~ r) ^ cos (|/3| 2 sin[2(A - r)<f>]) L n (4|/3| 2 sin 2 [(7V - r)<f>]) > 

is the normalization constant, and L n (z) is the Laguerre polynomial. An additional property of these superpositions 
was established in [11] for n = and 1 <C (e\(3\ 2 /N) N <C 4 W , where the interference effects approximately produce 
a particular Fock state; while for (e\/3\ 2 /N) N <C 1, an almost vacuum state is reached. To engineer (1) we initially 
prepare the trapped-ion state in |$(0)) = |n) ® | f) by means of the experimental techniques described in Refs. [2-5] 
(this procedure characterizes the first step of our scheme). With the help of the method established in [35] and used 
by Moya-Cessa et al. [10] for obtaining an arbitrary superposition of coherent states, the second step consists in the 
generation of superpositions of two displaced number states on a line. Now, considering the proposal of Duarte Jose 
and Mizrahi [11] for engineering circular states and adopting the motional state reached in the previous procedure 
as the initial motional state for this last step, we obtain, after a sequence of I cycles involving the application of 
laser pulses and no-fluorescence measurements, the state (1) in phase space. Furthermore, the total duration T of 
laser pulses employed in the sequence and the probability P|-(T) of getting the ion in the upper electronic state 
after I cycles are explicitly calculated. We also verify that the quantum interference effects between the TV = 2 i+1 
components of the motional state obtained in the third step decrease the values of Pj (T) when superpositions with 
N » 1 are regarded. Finally, assuming that the effective relaxation process of trapped ions can be described by the 
standard master equation for the damped harmonic oscillator [16, 25], we evaluate the time evolution of the Wigncr 
function associated to \^>^ (/?)). Following, this function is factorized into diagonal and nondiagonal terms which 
permits us to investigate, for example, the degradation of the undermentioned quantum interference effects through 
a quantitative measure of coherence introduced in [36] that characterizes the decoherence process of the nondiagonal 
elements of a density operator in the Fock-state basis. 

This paper is organized as follows. In Section II we adopt the method proposed by Wallentowitz and Vogel [35] 
in order to obtain a unitary time-evolution operator which permits one to produce superpositions of two displaced 
number states on a line. To engineer (1) we consider in Section III the proposal of Duarte Jose and Mizrahi [11] for 
obtaining circular states, and also determine the total duration of laser pulses and the probability of getting the ion 
in the upper electronic state during the construction process. In Section IV we employ the Weyl-Wigner formalism to 
investigate the degradation of the quantum interference effects between the 2 l+1 components of the motional state (1) 
through an effective relaxation process of the trapped ion. Section V contains our summary and conclusions. Finally, 
Appendix A describes the calculational details on the measure of coherence used in Section IV. 



II. ENGINEERING SUPERPOSITIONS OF TWO DISPLACED NUMBER STATES ON A LINE 

Let us consider a weak electronic transition of an ion which is bichromatically irradiated by two laser fields detuned 
to the first lower and first upper vibrational sidebands of the transition, respectively, with equal intensities. In the 
resolved-sideband and Lamb-Dickc regimes, the interaction Hamiltonian for the laser-assisted vibronic coupling can 
be written, in the interaction picture, as [35] 

H int = V/2 SI (<x_ + <t+ e-*) X e , (2) 

where il = 77A is the effective Rabi frequency on the first vibrational sideband with coupling constant A and Lamb- 
Dicke parameter r\. The electronic flip operators cr± and <x z describe the electronic transitions | J.) ^ | |) and satisfy 
the commutation relations [er + ,cr_] = er z and [cr z ,er±] = ±2cr±. The phase-rotated quadrature operator 

X e = a 6 + * 6 — = Qcosfl-Psintf {-tt < < tt) (3) 
v2 

represents the generalized CM position of the ion, being a* (a) the creation (annihilation) operator of vibrational 
quanta. Here, the dimensionless quadrature operators Q (positionlike) and P (momentumlike) obey the Weyl- 
Heisenberg commutation relation [Q,P] = il (for simplicity, we will fix ft = 1 throughout this paper). Furthermore, 
the phases p — \{pb + <fii) and 6 — -^(fb — <Pi) contain the phases ip^ and ip r of the lasers detuned to the blue (b) 
and red (r) sides of the electronic transition. In particular, when <p r = ip^ or <p T = p^, + ir, we obtain the operators 
X = Q and X_ 7r/2 = P. 

Using the interaction Hamiltonian (2), we can express the unitary time-evolution operator U(t) = exp (— itH int ) in 
a compact form as follows [10]: 

U(t) = cos (V2SUX0) - i (<r_ e iv + er + e"^) sin (V2OtX ) . (4) 
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FIG. 1: Plot of P T (/3) = |[1 +exp(~2|/3| 2 )L n (4|/3| 2 )] versus \/3\ £ [0,4] and different values of the excitation degree n, where 
the dot-dashed, dashed and solid lines correspond to n = 0, I and 2, respectively. 



This result allows one to determine the density operator p(t) = U(t)p(0)W (t) with p(0) = p v (0) <£> | |)(T I (being 
p v (0) the density operator for the CM motional state at time t = 0), and to prepare a superposition of displaced 
number states on a line in phase space. In fact, we are interested in generating superpositions of two displaced 
number states on a line using the present approach. For this purpose, the evolution operator (4) is applied on the 
state |$(0)) = |n) <g> | f), giving the following result: 

|$(/?)> - \ (|n, (5) + \n, -/3» ® | T> - \ e iv (|n, (5) - \n, -/?)) ® | |> , (5) 

where |n, ±/3) = D(±/3)|n) correspond to the displaced number states [37] and whose statistical properties were 
studied in detail by de Oliveira et al. [38], D(/3) = cxp (/3a* — (3*a) is the displacement operator, and /3 — iQte~ l9 . 
The procedure of measurement of the motional state was established in [2-5] , and it consists of collecting the emitted 
resonance fluorescence signal from the transition \d) «-> | j) (being \d) an auxiliary electronic state with width T on the 
order of r/27r « 20 MHz) by means of a laser strongly coupled to the electronic ground state during a specific period 
of time t. Following, we consider only those events where no fluorescence have been observed, since any spontaneously 
emitted photon will disturb the motional quantum state via recoil effects. At this point, it is important mentioning 
that the efficiency in collecting the fluorescence of the trapped ion is of order of 10~ 4 (i.e., about 10 4 photons have 
to be scattered by the cycling electronic transition to be detected). Thus, the time needed to detect the electronic 
state of the ion is approximately equal to Td ~ 200/xs [39, 40]. Now, if one considers the measurement time (rd) in the 
evaluation of the total time (r t ) necessary to prepare the target state, one obtains r t ~ r + 200^s. Consequently, the 
quantum state (5) is projected onto the excited state | f) and the resulting conditioned vibronic quantum state reads 

|$( / 3)) = 7 -^-(|n,/3) + |n,-/3))®|T), (6) 

being A/"i 2 ' the normalization factor. In particular, the probability P|(/3) for the occurrence of the no-fluorescence 

event is connected with the normalization factor through the relation |A/"i 2 ^| 2 Pf (/3) = 1/2, and its maximum and 
minimum points depend on the excitation degree n of the motional state described by Eq. (6). Figure 1 shows the 
plot of P^((3) versus \(3\ for n = (dot-dashed line), 1 (dashed line) and 2 (solid line), where we observe that the 
maximum of this function reaches P-\ ~ 0.68 when n = 2 and \(3\ ~ 1.27 (for instance, if one considers r\ w 0.1 and 
X/2tt s=a 1 MHz, this value of corresponds to T t « 202 /is which is greater than 27r/r w 0.05 fis). 



III. ENGINEERING SUPERPOSITIONS OF DISPLACED NUMBER STATES ON A CIRCLE 

In this third step, we adopt the procedure established by Duarte Jose et al. [11, 16] which is based on a Kerr-type 
interaction obtained through the interaction between the trapped ion and one pair of laser beams tuned in resonance 
with the electronic transition frequency In the Lamb-Dicke regime (for more details, see Refs. [41-43]), the carrier 
Hamiltonian of this system can be approximated as follows [11]: 



7i = Acr x cxp ( — - 



K 4 K 6 

1 - K 2 n H (n 2 - n) (n 3 - 3n 2 + 2n) + 

4 36 



Acr x (l - K 2 n) , 



(7) 
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where A is the effective Rabi frequency, n is the Lamb-Dicke parameter, n = a^a is the phonon-number operator, and 
1 is the identity operator. It is important mentioning that the validity of the Hamiltonian Ti. essentially depends on 
the condition (n 2 / A) (a^ 2 a?) <C (a^a), i.e., this condition must be satisfied in order to guarantee the validity of the 
approximation employed in (7). In addition, we also consider the vibronic quantum state (6) as an initial state for 
the third step of the protocol. 

To construct superpositions of displaced number states on a circle we have used the sequence outlined in Ref. 
[11] which consists of £ cycles involving the application of laser pulses and no-fluorescence measurements. In fact, 
each cycle consists in the application of one laser pulse with specific duration tk (in order to generate the required 
superposition) followed by one no-fluorescence measurement (this event assures no-recoil effects of the vibrational 
motion of a trapped ion, and maximizes the probability to realize successfully the target state). Now, if in a particular 
cycle of measurements a fluorescence emission is detected, the sequence must be stopped and repeated again. Thus, 
after I cycles of successfull measurements, the resulting conditioned vibronic quantum state becomes 



|*(ti + ■■■+*/)) =y/2M%> 



n<T |U(t*)| T) 



L fc=i 



{\n,0) + \n,-0))®\ T) 



where 



U(i fc ) = cxp [-itk (Al - An) a x ] (A = k 2 A) 



(8) 



(9) 



is the unitary time-evolution operator at time tk associated to the carrier Hamiltonian 7i. To engineer the ion CM 
motional state as Eq. (1), we need to adjust the phases in (8) of the displaced number states putting them evenly 
distributed around the circle: this fact is only possible when the duration of the fcth evolution pulse is given by 
tk = ^ I (2 fe+1 A). Furthermore, we have chosen conveniently the Lamb-Dicke parameter as n 2 = (n + 2 e+2 ) , which 
permits us to write the superposition (8) in the simplified form 

|*(ti + --- + t/)> = |*^ +1) ()9)>®|T> (10) 

since <j> = n/2 e+1 , 6 r = 2i:r/2 t+1 with r = 1, . . . , 2 e+1 , and N = 2 e+1 are fixed a priori. The final adjustment of 
phases involved in this process is reached when the phases ip T and ipt, of the lasers described in the previous section 
satisfy the relation ip T = <^ b + 7t/4. Consequently, the normalization constants M$ and N^ 2 - 1 can be connected by 

means of the equality M ( n ] = 2( /+1 >/ 2 N< l 2<+1 >. 

The total duration T = t\ + ■ • ■ + 1( of laser pulses employed in the sequence, the total time T t necessary to prepare 
the target states (10), and the probability Pj(T") of getting the ion in the upper electronic state during the £ cycles, 
i.e., 



7T 



2A 

T t - T + £r d 



(n + 2 e + 2 ) (1 - 2- e ) 



Tt 



(11) 

(12) 



and 



P T (T) 



1 



1 



2^+i 1 2 2e+1 



r exp 



l/3| 2 sin 2 (^ 



cos 



' sin 



2nr 
2l+T 



4|/3| 2 sin 2 (, 



(13) 



permit us to characterize completely the construction process. It is important mentioning that N^ 2 + •* and P-\{T) are 

connected through the relation |N^ 2 ^\ 2 P^(T) = 2~ 2 ( £+1 ). As an application of the results obtained until the present 
moment we consider the engineering of superpositions with = 4 and 8 displaced number states, which correspond to 
different sequences involving I = 1 and 2 cycles each one, and analyze their respective success through the probability 
(13). Figure 2 shows the plot of P^(T) versus \/3\ for (a) £ = 1 and (b) £ = 2, with n = (dot-dashed line), 1 (dashed 
line) and 2 (solid line) fixed in both situations. In Fig. 2(a) the maximum of this function reaches Pj w 0.37 (0.34) 
when n = 1 (2) and \/3\ w 1.65 (1.28), with k w 0.33 (0.32) and T t w 401.93 /j,s (401.85 /j,s) if one considers 2tt/A w 1 
/is. On the other hand, in Fig. 2(b) this maximum reaches Pj w 0.25 (0.20) for n = 1 (2) and \f3\ w 1.96 (3.03) with 
k w 0.24 (0.23) and T t w 604.5 ^s (605.5 /zs). Furthermore, note that for n = 2 and \(3\ w 1.27 the probability P T (T) 
is approximately equal to 0.34 (0.18) when the superposition has N = 4 (8) states. In fact, this value depends on 
the number of cycles involved in the sequence and decreases when superpositions with N 3> 1 are regarded (compare 
these values with that obtained in Fig. 1). The explanation of this result is associated with the quantum interference 

effect between the 2 e+1 components of the motional state i- c -, high values of N lead us to obtain a large 

number of components interfering with each other, and this interference decreases the value of Pf(T). 
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FIG. 2: Plot of -Pf(T) versus \/3\ £ [0,4] for superpositions involving (a) N = 4 (one cycle) and (b) N = 8 (two cycles) 
displaced number states. The total duration of laser pulses employed in both sequences for each excitation degree n is much 
more than that proposed by Ref. [44], where the number states were generated through a sequence of quantum nondemolition 
measurements on a thermal and a coherent initial state. Here, the dot-dashed, dashed and solid lines correspond to n = 0, 1 
and 2, respectively. 



IV. THE DEGRADATION OF THE QUANTUM INTERFERENCE EFFECTS VIA WIGNER 

FUNCTION 



In this section, we adopt the procedure established by Refs. [16, 42] and assume that the effective relaxation process 
of trapped ions can be described in the framework of the standard master equation for the damped harmonic oscillator 

[45], 



dPv 
dt 



= -iw [a f a, p v ] + 7 (n + 1) (2ap t ,a t - a<!ip v - p v a<a) + jn (2aV„a - aa f p v - p,„aa f ) , 



(14) 



where a (a*) is the annihilation (creation) operator associated with the oscillatory motion of frequency wo in a 
one-dimensional harmonic trap, p v describes the density operator for the CM motional state at time t, n is the 
equilibrium mean number of motional quanta in the reservoir, and 7 is a positive relaxation rate of the energy to 
thermal equilibrium (the connection between master equation and averaged interferometer approach in the trapped 
ion context was addressed in Ref. [25]). If one considers the Wigner representation, this equation is equivalent to the 
Fokker-Planck equation, 



dW 



d d 

— (79 - u p) + q^(1P + ^oq) + l(n + 1/2) 





d 2 V 


\8q 2 


F dp 2 ) _ 



W, (15) 

for the time-dependent Wigner function W(p,q;t) 1 whose solution can be written as an integral equation [36] 

(16) 



/OC 
K(p, q; t\p>, q'; 0)W(p', q'; 0) dV (dV = dp' dq' 
-00 

with the kernel 

K(p, q; t\p', q 1 - 0) - [tt(1 + 2n)u]- 1 cxp {- [(1 + 2n) U f 1 [( Pt - e^p') 2 + (q t - 



depending on the time variable and reservoir parameters, p t = pcos(Ld a t) + qsin(ui a t), q t — qcos(uj t) — psin(uj t), 
and u{t) = 1 — e~ 27 * (this function was denominated as 'compact time' in Refs. [16, 36]). 

(2 e+1 ) 

In order to calculate the time evolution of the Wigner function associated to the motional state (/?)), firstly 

we substitute into the integrand of Eq. (16) the initial Wigner function [34] 



W n (p',q';0) = ^|Nf 



J )|2 



^ cxp (— d\ rr ) L n (2$R rr ) 



r=l 



-1 2* 



+ 2 ^ ^2 exp[-Re(m rs )]coa[lm(m rs )}L n [2Re(m rs )} 



(17) 



= 1 r=s+l 
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with 



*r.(0) = [(<?' + W') - V2 Pr] [(</ + ip') - V2 ft] * + \f3\ 2 - 0rfi 
Then, carrying out the integrations in the variables p' and q', we get 



W n (p, q;t) 



(-1)" |N<? 



2^)12 



7r 1 + 2nu 



-1 2 e 



1- 2(n+ l)u 



1 + 2nw 



n I 2 



^ exp (-$rr) L n (20 rr ) 



+ 2 Re[exp(-# rs )L n (20 rs )] 

s— 1 r— s+1 



(18) 



where 



Re [&•.(*)] 
Im [£„(*)] 
Re[C5 rs (i)] 
Im [«„(*)] 



Re[fH rs (t)] 2[(l + 2n) 



1 + 2nu 



+ 



1 + 2nu 



■ sin 



7r(r — s) 



Im[9l rs (t)] [(l + 2n) U ]|/3| 2 . 

■H : : Sin 



1 + 2nw ' 1 + 2nu 
(l-u)Re [m rs (t)} 



2 e+i 
2ir(r — s) 



2 e+i 

2[(l + 2n)w] 21 " 2 



sin 



(1 + 2nu)[l - 2(n + l)u] (1 + 2nu)[l - 2(n + l)u] 
[(1 + 2fi)u]Im [m rs {t)} (l-u)[(l + 2n)u]\l3\ 2 ^ 



ir(r — s) 



2 e+i 

2ir{r — s) 
2 l + l 



and 



(l + 2nu)[l-2(n+l)u] (1 + 2rat)[l - 2(n + 
9tr.(*) = [(Qt + ipt) - a/2(1-u) /3 r l f(q t + ipt) - V2(l-u) /?J * + (|/?| 2 - /?r/3 s *) (1 - «) • 



Note that W n (p, q; t) is factorized into diagonal and nondiagonal terms. This permits us, in particular, to investigate 
the degradation of the quantum interference effects among the 2 e+1 components of the motional state represented 
by the initial Wigner function W n (p,q;0)- Similarly, Chountasis and Vourdas [46, 47] have employed the same 
factorization for the Weyl and Wigner functions associated with a superposition of m quantum states |sj), and 
showed that the nondiagonal terms describe the interference effects between the states |sj). To illustrate these results, 
in Fig. 3(a) we have plotted the three-dimensional picture of W n (p,q;0) versus p and q (containing both diagonal 
and nondiagonal terms) for n — 2, \/3\ — 3.03, and 1 — 2 fixed; while (c) and (e) correspond to diagonal and 
nondiagonal terms, respectively. The influence of these terms on the shape of the Wigner function (17) leads us to 
confirm the results previously obtained by Chountasis and Vourdas since the nondiagonal term (c) is responsible for 
the interference pattern observed in (a). Now, the degradation of this pattern for t > is connected with the effective 
relaxation process under consideration. To illustrate this point, we have plotted in Fig. 3(b) the three-dimensional 
picture of W n (p, q;t) for the same parameter set used in the previous figure, with addition ofwo/7 — 1, n = 1 and 
7< = 0.1. Furthermore, Figs. 3(d) and (f) represent the diagonal and nondiagonal terms, respectively. From the 
comparison between Figs. 3(a) and (b) we can perceive that the quantum inteference pattern present in the first 
picture has disappeared in (b), and this fact is associated with the decoherence effect on the nondiagonal elements 
of the density operator p v it) (here mapped into the nondiagonal term of W n (p, q; t) and pictured through figure 
3(f)). It is important mentioning that (d) also has a significant contribution to the shape of (b), and when w 1 
this contribution is dominant if one compares with that obtained from the nondiagonal term. Similar results can be 
reached if — 0.1 and n » 1, since the equilibrium mean number of quanta in the reservoir represents a scale factor 
for the compact time u(t). 

The measure of coherence C(t) was introduced in [36] as a quantitative measurement which characterizes the rate of 
decoherence in the Fock-state basis. Since decoherence can be interpreted as the disappearence, with time progression, 
of the nondiagonal elements associated to the density operator p v (t) , we will use this measure in order to explain the 
pattern observed in the diagonal and nondiagonal terms of W n (p, q;t). For this purpose, we define the normalized 
measure of coherence through the expression (for more details, see Appendix A) 



<#>(t) 



m1 £) (o)-a^(o) 



(19) 
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FIG. 3: The pictures (a,c,e) represent the three-dimensional plot of W n {p, q\ 0) versus p and q for n = 2, \(3\ — 3.03, and 1 = 2; 
while (b,d,f) correspond to W n (p, q; i) with W0/7 = 1, n = 1, and 7t = 0.1. The interference pattern observed in (a) (containing 
both diagonal and nondiagonal terms) is a direct consequence of the nondiagonal term (e), and its disappearance in (b) is 
connected with the effective relaxation process under consideration. 



where the 'total purity' fin (t) and 'diagonal purity' X n (t) can be written in terms of the initial Wigner function as 
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FIG. 4: Plot of Cn\t) versus u(i) for \/3\ — 3.03, 1 = 2, and n = 1 fixed, where the dot-dashed, dashed, and solid lines 
correspond to n — 0, 1, and 2, respectively. Note that in a short period of time (-yt G [0,0.1116]) and high values of £ and n, 
the measure of coherence suddenly goes to zero. 

follows: 

/oo poo 
/ E(p', q ',p",q";t)W n ( P ', q ';0)W n ( P ",q";0)dr'dT" , (20) 
- oo J — oo 

/OO /" oo 
/ A(y, t? ',p", t7 ";t)T47 n (p', (7 ';0)T47 n (p", t? ";0)dr'dr", (21) 
-oo J — oo 

being S and A given by Eqs. (A7) and (A10), respectively. The functions fin\t) and Xn\t) were explicitly calculated 
in Appendix A, and their analytical results for n = corroborate that obtained by Souza Silva et al. [16]. Figure 4 

represents the plot of C„\t) versus the compact time u(t) = 1 — e~ 2lt for |/3| = 3.03, £ = 2, and n = 1 fixed, where 
the dot-dashed, dashed, and solid lines correspond to n = 0, 1, and 2, respectively. Note that for u = 0.2 we obtain 
ef ] w 0.0443, Cf ) w 0.0438, and d, 2) w 0.0076; while for n = 2, the measure of coherence goes to cf ] w 0.0428, 

~ 0.0247, and ~ 0.0076 at the same compact time. Otherwise, if one considers u = 1 in both situations, dp 
goes to zero for any values of £ and n. Thus, high values of the parameters i and n considerably decrease (increase) 
the measure of coherence (decoherence process) in a short period of time and this fact could explain the disappearance 
of the quantum interference pattern observed in Fig. 3 [48] . 




V. SUMMARY AND CONCLUSIONS 



In this paper, we have combined different theoretical approaches in order to engineer superpositions of displaced 
number states on a circle in phase space as target states for the center-of-mass motion of a trapped ion. The total 
duration T of laser pulses employed in the process and the probability Pj (T) of getting the ion in the upper electronic 
state were explicitly calculated and analyzed. In particular, we have verified that (i) the quantum interference 
effects among the N = 2 l+1 components of the motional state described by Eq. (1) decrease the values of P-\(T) 
when superpositions with N 3> 1 are regarded, and (ii) the Lamb-Dicke parameter k essentially depends on the 
number £ of cycles involved in the sequence and the excitation degree n of the motional state. Furthermore, we 
have also investigated the degradation of the quantum interference effects via the Wigner function and showed that 
(iii) these effects basically depend on the nondiagonal term of the Wigner function at time t — 0, (iv) the quantum 
interference pattern present in W n (p, q; 0) disappears for t > and this fact is associated with the decoherence process 
on the nondiagonal elements of the density operator p v (t), and (v) high values of the parameters £ and n increase 
the decoherence process in a short period of time. Summarizing, the work reported here is clearly the product of 
considerable effort and constitutes an marginally original contribution to the wider field of quantum state engineering. 

Recently, Lvovsky and Babichev [49] have synthesized the displaced Fock states of the electromagnetic field by 
overlapping the pulsed optical single-photon Fock state with coherent states on a high-reflection beam splitter and 
showed its nonclassical properties (such as negativity of the Wigner function and photon number oscillations) through 
a complete tomographic reconstruction. However, the nonunitary quantum efficiency of the homodyne detector, the 
dark counts of the single-photon detector, and the impurity of the optical mode of the conditionally prepared photon 
represent important restrictions on the preparation and measurement of the Fock state. In this sense, experiments 
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involving trapped ions are a great laboratory in the construction process of nonclassical states since the decoherence 
time is the longest if one compares with that obtained from experiments for trapped and travelling nonclassical 
states of the electromagnetic field. In conclusion, we believe that the results obtained in this paper can motivate 
the generation of new nonclassical states in future experiments on trapped ions and to contribute significantly to the 
study of quantum interference effects in different physical contexts. 
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APPENDIX A: THE MEASURE OF COHERENCE 



Dodonov et al. [36] have introduced two quantitative measures which characterize the rates of decoherence and 
thermalization of quantum systems, and studied the time evolution of these measures in the case of a quantum 
harmonic oscillator whose relaxation process is described in the framework of the standard master equation. In 
particular, the measure of coherence C(t) was defined by the authors through the expression 

C[t) ~ M0) - A(0) ' (A1) 

where the functions u.(t) = Tr [p 2 (t)] and A(t) = Tr [p^(i)] (Pd(t) — J2nen Pn(t)\n) (n\ correspond to the diagonal 
part of the density operator p(t), being P n (t) = (n\p(t)\n)) connected with Wigner function W(p,q;t) and phonon 
(photon) distribution function P n (t), as follow: 

/oo 
[W(p,q;t)] 2 dT (dT = dpdq) (A2) 
-OO 

and 

oo 

\(t) = J2l p n(t)} 2 • (A3) 

n=0 

In addition, the phonon (photon) distribution function can also be obtained by means of the auxiliary relation 

Pn(t) = / W n (p, q)W(p, q; t) dT , (A4) 

J — OO 

being 

W n (p, q) = 2(-l)" cxp [- (p 2 + q 2 )] L n [2 (p 2 + q 2 )] 

the Wigner function associated to the number state at time t = 0. Thus, if one knows the Wigner function associated to 
the density operator p(t), then the functions 'total purity' /i(t) and 'diagonal purity' X(t) can be promptly calculated. 
Note that C(0) = 1, and C(t) =0 for any completely incoherent state without nondiagonal matrix elements in the 
energy basis (provided that initially at least one nondiagonal element was different from zero). 

Now, substituting the solution (16) of the Fokker-Planck equation for the harmonic oscillator into Eqs. (A2) and 
(A4), we obtain 

/OO /'OO 
/ Z(p>, q',p", q"; t)W(p', q'; 0)W(p", q"; 0) dr'dr" , (A5) 
-oo J — oo 

/oo 
JC n {p',q';t)W(p',q';0)dT' , (A6) 
-oo 



11 



with 



Z(p',q',p",q";t) = [(1 + 2n)u] 1 exp 



(i- M ) [(p' -p") 2 + W -q"Y] 

2(l + 2n)u 



(A7) 



and 



IC n (p',q';t) 



2(-iy 



1 + (1 + 2n)u 



1 - (1 + 2n)u 



1 + (1 + 2n)u 



exp 



(l-u)(p' 2 + q' 2 ) 
l + (l + 2n)u 



x L r , 



2(1 - u) (p' 2 + q' 2 ) 



[l-(l + 2n)u][l + (l + 2n)ti] 
Consequently, the function A(t) can also be determined through the equation 

/oo poo 
/ A(p', q',p", q"; t)W(p', q'; 0)W(p", q"; 0) dr'dr" , 
-oo J — oo 

where A(p' , q',p", q"; t) is given by 

(1 - u)(p' 2 + q' 2 + p" 2 + q" 2 ) 



(A8) 



(A9) 



A(p',q',p",q";t) = [(1 + 2n)u]- 1 exp 
x I a 



2(1 + 2n)u 

(1-u) [(p> 2 + q > 2 )(p" 2 + q " 2 )] 1/2 
(1 + 2n)u 



(A10) 



being I v (z) the modified Bessel function of the first kind [50]. This expression was obtained with the help of the 
intermediate relation [51] 

n— ^ ' 

with \z\ < 1 and v > — 1. Since the initial Wigner function does not depend on the undermentioned relaxation process, 
Eqs. (A5) and (A9) represent an alternative way to the calculation of the total and diagonal purities. 

In many situations of practical interest both purities can be established rather easily. For instance, if one considers 
the initial Wigner function (17) in Eqs. (A5) and (A9), we get 



z(x + y) 




'2{xyzfl 2 ' 


1 - z 




l-z 



2iN<rv a 



a 



In 



(xy) n 



dx n dy n xA y + (1 — u)y 



R(x,y;t) 



x=y=l 



and 



WW 



2|Nf 



')|2 



nl 



Q2n 



(xy) n 



dx n dy n A+B+ - A-B. 



3(x,y;t) 



(AH) 



(A12) 



x = y = l 



where 



2 f+i 2 e+i 

R(x,y;t) = ]T cxp[-2(^-£-)|/3| 2 ] , 

p,q—l r,s—l 

r B+A rs V rs + A+B pq V pq , 2(1 - u) (B+B-XJ 2 rs + A+A-V. 



3(x,y;t) = 



S H exp - 

p,q—l r,s—l 



A+B, 



+ 



A+B+ (A+B+-A-B-) 



X In 



8(1 -u)U ra V p ,|/ 
A+B+ - A-B- 
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with 



€£(x,y;t) = [xk y + (l-u)y] 1 |a„U^ + 2(1 - u) cos 



7r(r + s — p — q) 
W 1 



U rs V P9 + A x V 2 pq , 



D r pq (x,y;t) = cos 

\j rs (x) = cos 

Vp g (y) = cos 

A rs (t) = cos 

Bpg(i) = cos 



7r(r — s) 
7r(r — s) 
7r(p - gr) 



U rs + cos 



2 i+i 
n(r — s) 

tt( P - q) 



x + i sin 



y + i sin 



' n(p-q) 
2 £+l 

n(r — s) 



2 i+i 

tt(p - q) 

2 £+l 



2 e+i 



(1 — u) — i sin 
(1 — it) — isin 



7r(r — s) 
q) 



2 i+i 



[1 + (1 + 2n)u] , 
[1 + (1 + 2n)u] , 



>l±(ar;i) = [1 ± (1 + 2n)u]ar ± (1 - u) , 

B±(y;t) - [l±(l + 2n)u]y±(l-u), 

A x (f) = 2(1 + 2n)ux + (1 - u) , 

Aj,(t) = 2(1 + 2n)uy + (1 - u) . 

Furthermore, the phonon distribution function is given by 



... (_-\\n+m fin Am n 

pWm = 2iN (2 + )i 2 L_y • A - a: i 



a;=l 



(A13) 



where 



I m (a;;f) = ex P 



r,s— 1 



2A rs U rs |/3| 2 



4(1 -u)XJ 



^l/?| 2 



A+A- 



It is important mentioning that Eqs. (All)-(A13) were calculated by means of the parametric representation for the 
associated Laguerre polynomial [52], i.e., 



x n+a e~ xz 



The total and diagonal purities determined in this appendix corroborate that obtained by Souza Silva et al. [16] for 
n = 0. 
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